# -*- coding: utf-8 -*-
"""
Created on Sat Sep  7 20:23:32 2024

@author: LENOVO
"""

import matplotlib.pyplot as plt
import numpy as np
from sympy import *
from scipy.optimize import root, fsolve
import pandas as pd
import math

H=1.7
a=16*H  #m
b=H/(2*np.pi) #m
vh=1 #m/s
L0=(341-27.5*2)/100
Lb=(220-27.5*2)/100

def rt(x):
    if x<=32*3.14159265358:
        r=((a-b*(x)))
    else:
        r=1
    return r

R=[]
for i in np.linspace(0, 64*np.pi,10000):
    r=rt(i)
    R.append(r)
theta=np.linspace(0, 64*np.pi,10000)



ax=plt.subplot(111, polar=True)
ax.set_theta_direction(-1)
plt.plot(theta,R,lw=1,c='r', label='各节点')
plt.legend()
plt.title('全过程')
plt.grid()  
plt.show()